GENERICITY OF BLACKHOLE FORMATION IN THE GRAVITATIONAL COLLAPSE OF 
HOMOGENEOUS SELF-INTERACTING SCALAR FIELDS 
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ABSTRACT. The gravitational collapse of a wide class of self-interacting homogeneous scalar fields models 
is analyzed. The class is characterized by certain general conditions on the scalar field potential, which, in 
particular, include both asymptotically polynomial and exponential behaviors. Within this class, we show 
that the generic evolution is always divergent in a finite time, and then make use of this result to construct 
radiating star models of the Vaidya type. It turns out that blackholes are generically formed in such models. 



1. Introduction 

Scalar fields attracted a great deal of attention in cosmology and in relativistic astrophysics in the last 
twenty years. It is, indeed, worth mentioning the fundamental role that they play in modeling the early 
universe and, on the other hand, the relevance of the self-gravitating scalar field model as a test-bed for 
the Cosmic Censorship conjecture, that is, whether the singularities which form at the end state of the 
gravitational collapse are always covered by an event horizon, or not || |l7| . In particular, in an important 
series of papers, Christodolou studied massless scalar field sources minimally coupled to gravity. For this 
very special case the formation of some naked singularities has been shown, but it has been also shown 
that such singularities are non generic with respect to the choice of the initial data f|]. 

The massless scalar field is a useful model in studies of gravitational collapse since its evolution 
equation in absence of gravity - the massless Klein-Gordon equation on Minkowski spacetime - is free of 
singular solutions. However, of course, the general case in which the field has non-zero self-interaction 
potential (which of course include the massive scalar field as a quadratic potential) is worth inspecting 
from the theoretical point of view. Further to this, in recent years, interest in models with non vanishing 
potentials aroused from string theory, since, for instance, dimensional reduction of fundamental theories 
to four dimensions typically gives rise to self-interacting scalar fields with exponential potentials, coupled 
to four-dimensional gravity. In this context, the case of potentials which have a negative lower bound 
plays a special role, since such potentials generate the anti De Sitter space as equilibrium solutions 

Bill- 

As a consequence of their "fundamental" origin, exponential potentials have been widely considered 
in the cosmologically-oriented literature (see[p5| and references therein) also with the aim of uncovering 
possible large-scale observable effects |20|, |2 jj, p6| ], A few examples of non trivial potentials (to be 
discussed in full details in the sequel) have been studied in the gravitational collapse scenario as well. In 
these papers formation of naked singularities has been found [T|]. 

So motivated, in the present paper we study the dynamical behavior of FRW (spatially flat) scalar 
field models, imposing only very general conditions on the scalar field potential; our hypotheses essen- 
tially reduce to ask the potential to be bounded from below and divergent when the field diverges. The 
evolution of collapsing solutions or, equivalently, time backwards expanding solutions, has been studied 
previously in the literature mainly making use of dynamical system approaches. From this point of view 
Foster's work is worth noticing; indeed, in [Q] existence of past singularities of expanding, spatially flat 
cosmologies are found in terms of sufficient conditions on the (nonnegative) potential V(<f>) - in partic- 
ular, the proof relies on the existence of a coordinate transformation that brings [<f>o, +00 [ in ]0, e] , and 
satisfies some other regularity conditions together with the transformed potential (see [0, Definition 2]). 
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In these cases, the behavior is proved to be essentially similar to the behavior in the massless case, and 
so a past singularity exists and it is stable with respect to the choice of the data. 

Other worth-mentioning works are those related to the case of non-minimal coupling, such as 
in which flat FRW model are studied using conformal transformations that bring the models into the 
standard form, and [|[, where singularity formation for non minimal coupling models is analyzed for 
some specific form of the potential. 

As mentioned above, in this paper we consider a class of potentials £ which is quite more general 
with respect to those previously considered. In particular, we weaken assumptions about the behavior of 
the potential at infinity (which are necessary in the dynamical systems approach) thanks to an accurate 
use of the Implicit Function Theorem. Indeed a suitable estimate of the involved neighborhoods is used 
to obtain a result of existence, uniqueness, and instability of solution for a particular kind of ordinary 
differential equation of first order to which the standard theory does not apply straightforwardly (Theorem 
3.6). Another key difference with previously published results is that the potentials may assume negative 
values. Non negativity of the potential means that the total energy of the scalar field has a positive 
lower bound (using this constraint many results for scalar field cosmologies have been actually found 
[0, [l2| £4]]; see also Jll|]); however, in many issues and especially when string theory comes 

into play, it becomes relevant to inspect spacetimes in which the potential still has a lower bound but 
this bound is negative; under certain conditions indeed, potentials of this kind generate solutions with 
positive total energy Jl6|], Since, in particular , a constant negative potential generates an "equilibrium" 
state which is just the Anti de Sitter solution, within this framework models which are asymptotically 
anti-De Sitter can be inspected [£l| ^ . 

Within the class £ of potentials (to be defined in precise technical terms in the section below) we 
prove that - except for a zero measure set of initial data - all the solutions become singular in a finite 
time. Considering such solutions as sources for models of collapsing objects (with generalized Vaidya 
solutions as the exterior) we study the nature - blackhole or naked singularity - of such objects. It turns 
out that, except possibly for a zero-measure set of initial data, all these "scalar field stars" collapse to 
black holes. It follows that the existing examples of naked singularities which occur in the literature must 
be non-generic within the set of possible solutions, and indeed we show in details how this occurs. 

Homogeneous scalar field models with self-interacting potentials thus respect the weak cosmic cen- 
sorship hypothesis, at least for a wide class of potentials. Of course, however, the problem of extending 
this result to the non-homogeneous case remains open. 



2. Collapsing models 

We are going to consider here homogeneous, spatially flat spacetimes 

(2.1) g = -dt (g) d< + a 2 {t) [dx 1 ® dx 1 + dx 2 ® dx 2 + dx 3 da; 3 ] 

where gravity is coupled to a scalar field <fr, self-interacting with a scalar potential V((f>). Assigning this 
function is equivalent to assign the physical properties of the matter source, and thus the choice of the 
potential can be seen as a choice of the equation of state for the matter model. Therefore, as occurs 
in models with "ordinary" matter (e.g., perfect fluids), the choice of the "equation of state" must be 
restricted by physical considerations. Stability of energy, in particular, requires the potential to have a 
lower bound. We shall thus consider here potentials of class C 2 (M) bounded from below; the lower 
bound of the potential need not be positive however, so that class of solutions which include AdS minima 
can be considered. Further, we need a set of technical assumptions. 

We say that V belongs to the set £ if the following conditions are satisfied: 

Al) (Structure of the set of critical points) The critical points of V are isolated and they are either 
minimum points or non-degenerate maximum points 

A2) (Existence of a suitable bounded "sub-level" set) There exists a < (3 £ R such that 

<P>(3=> V'(4>) > 0, <j) < a => V'((f>) < 0. 
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Moreover, setting V* = max[ a ^j V, it is 

lim V(<j>) > K, lim V{<p) > V.. 

*— oo (p — >-\-oo 

A3) (Growth condition) The function defined in terms of V and its first derivative V' as 

{22) uW = ^mr 

satisfies to 

(2.3) limsup|u((/))| < 1, 

(2.4) 3 lim u'((f>)(=0). 

— >±00 

We stress that the set € defined in terms of the above conditions is sufficiently large to contain all 
the potentials which have been considered in physical applications so far. For instance, it contains the 
standard quartic potential (V(4>) = ^m 2 ^ 2 + A 2 ^ 4 ) and, more generally, all the potentials bounded from 
below whose asymptotic behavior is polynomial (i.e. <fi 2n for \(f>\ — > oo, n > 1). 

Once a potential has been chosen, the Einstein field equations allow, in principle, to obtain the dynam- 
ics of the coupled system composed by the field and the metric. These equations are given by 

3o 2 

(2.5a) (Gg = 8^T °) : ~ = -^ 2 + 2V(cf>)), 

a z 

(2.5b) {G\=^Tl): -^±^ = i< p 2 -2Vm, 

the Bianchi identity is of course, due to the Noether theorem, equivalent to the equation of motion for the 
field: 

(2.6) = -20 U + V'{4>) + Z-j^j = 0. 

Denoting the energy density of the scalar field by 

(2.7) e - 2 + 2V(<j>), 
We will actually consider the following system: 

2 

a \ e 



(2.8a) 

(2.8b) <t> + V\<f>) =-3 



a J 3 

a 
a 



Of course, it is mandatory to extract the square root in the first equation. Since we are considering the 
collapsing model we choose d to be negative and therefore our final system is composed by 

(2.9a) £ = -JI. 

a y 3 

(2.9b) 4> + V'(4>) = V3y^cj>. 

Let us observe that, using (2.7) and (2.9b), the following relation can be seen to hold: 

(2.10) e = 2\/3Ve0 2 , 

that is e(t) is monotone. Observe that the two cases - expansion and collapse - are connected by the 
fact that, in the expanding case, there might be the possibility of reaching a vanishing value of e in a 
finite time to. If this occurs, the model will be ruled from to onward by the equations for the collapsing 
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situation. Actually, it will be proved that all such solutions are divergent at some finite time for almost 
every choice of the initial data, thus yielding a singularity as well. 



In what follows, we are going to focus on solutions of the equations (2.8a)-(2.8b) - or, better, (2.9a)- 
( 2.9b| ). Actually, it can be proved [ |To| ] that, if a ^ for alH > and (a, (j>) are C 2 functions that solve 



(2.5a)-(2.5b) (with a non-constant 4>{t)) then they are solutions of (2.8a)-(2.8b) as well. To prove the 
converse, it suffices to derive (2.8a) w.r.t. time and use (2.8b) to deduce that if (in some interval [a, &]), a 
is not identically zero - that is, e ^ - then solutions of (2.8a)-(2.8b) are also solutions of (2.5a)-(2.5b). 



Finally, the unique case in which e(i) = 0, Vt > and the Einstein Field Equations are satisfied is given 
by the trivial solution a(t) = ao, <fr(t) = 4>o, with V(<f>o) = 0. 



3. Qualitative analysis of the collapsing models 



The aim of the present section is to study the qualitative behavior of the solutions of ( 2.9a )-( 2.9b ) 
whenever V <G (£. The initial data (on t = 0, say) for the system can be fully characterized in terms of 
the values of the scalar field <fio — 0(0) and of its first derivative. In what follows we will be mainly 
interested in constructing models of collapsing object and therefore we consider only solutions which 
satisfy the weak energy condition on the initial data surface. Therefore, we assume 



(3.1) 



e = e(0) = 2 (O) + 2V(4>(0)) > 0. 



Given a solution 4>(t), let I C [0, +oo) be the maximal right neighborhood of t = where the solution is 
defined, and call t s = sup I. The main result of this secti on i s that the solution cj>(f) diverges in a finite 
time for almost every choice of the initial data satisfying (3T ), that is t s is finite for almost all solutions. 
This result will in turn allow us (next section) to construct models of collapsing objects by matching these 
solutions with suitable Vaidya metrics. 

To proceed further, we start by considering only data such that eo is "sufficiently large", more precisely 



(3.2) 



e > 2V* 



where, for each choice of the potential V in £, the constant V * is the chosen as the same which appears in 
(A2) . Once the result will be proven for these data, we will extend its validity by removing the restriction 



We start our study with the following lemma. 
Lemma 3.1. If4>{t) is bounded, also <p(t) is bounded, andt s 



-co. 



Proof. Let K a bounded set such that <fi(t) e K,Vt E I, and argue by contradiction, assuming </> not 
bounded, and supposing that supj <p(t) — +oo - the same argument can be used if <f> is unbounded only 
below. Let t such that 

s u Pif V'(<j)) 



therefore (2.9b) implies that (j>(t) > 0, and this means that <j> is increasing in t, so ( 2.9bj ) shows that 
(j) > 0, Vt > i, and then <fi is eventually increasing, namely lim t _^ ts <fr(t) = +oo. This shows that 
t s < +oo, otherwise it would be <j)(t) — 0o = J Q 0(' r ) dr, which would diverge as t — * +oo. Then 
t s € R. 

With the position 



4>{t) 



equation (2.9b) implies 



A(t) 



vm)) 
m 2 




2V(4>(t)) 



V3, 



COLLAPSE OF HOMOGENEOUS SCALAR FIELDS 



5 



and then A(t), near t = t s , behaves like v3(*s — t). Since (f){t) positively diverges, lim t ^ ts (f)(t) 
exists, and is finite since <f>(t) is bounded by hypothesis. But the quantity 



- J 0(r)dr = J ^-ydr 

diverges, which is absurd. This shows that </>(t) must be bounded too. 

To prove that t s = +00, we proceed again by contradiction. Let tk — * t s be a sequence such that 
the sequence (<fi(tk), <j>(tk)) converges to a finite limit Solving the Cauchy problem (2.9b) 

with initial data (</>*, shows that the solution <f)(t) is C , the solution could be prolonged on a right 
neighborhood of t s . □ 

Proposition 3.2. The function <fi(t) is unbounded. 



Proof. Let by contradiction <f>(t) be bounded. Then, by Lemma 3.1, \4>{t)\ < M for a suitable constant 
M, and i s = supl = +00. In particular V(cf>(t)), V'((f>(t)), e(t),cf>(t) are bounded also. Then Lemma 
A.l in Appendix ^jsays that lim t ^ +00 4>{t) = 0, that implies lim t ^ +oc 2V((f>(t)) = lim t ^ +co e(t) > 
eq > 2V*, and so V(4>(t)) > V* for any t sufficiently large. Then, by assumption A2), for large t, <j>(t) 
moves outside the interval [a, (3] described in assumption A2), namely in a region where V is inve rtible . 
Since V(<j)(t)) converges this happen also for cj>{t) which converges to a point 0* $]a, [3[. Then by ( 2.9b ), 
ef>(t) -> -V'(0*) ^ (since -> 0), that is absurd. Then limsup t ^ ts \cj)(t)\ = +00. □ 

We have shown that (j>(t) is not bounded. Now we want to show that, actually, </>(i) monotonically 
diverges in the approach to t s . 

Remark 3.3. We observe that the quantity 

(3.3) " pit) := 2 -mi, 
satisfies the equation 

(3.4) p = 2 py/{l + p) (u(<P) y/T+p - sgn (0)) . 

Proposition 3.4. Let V G <£, and let <j> : I C R — » 1R f/ze solution of ( {2.9b ). </>(£) i's eventually 

nonzero, and 

lim |0(t)| = +00. 

£ — *sup I 

Prao/ Note that 4>{t) = imply 7(0(t)) = §e(i) > ±e > V*, so (f> g [a, /?]. Moreover if (j>{t) > (3 , 
it is V'((p(t)) > so 4>(t) < 0, while if cf>(t) < a it is ip(t) > 0, Therefore if by contradiction, <f>(t) = 
for infinitely many values of t, we see that cj>it) crosses the interval [a, (3} infinitely many times. 
Then there exist sequences o~k < Tk going to t s such that 

(f>(a k ) = a, <j>{r k ) = f3, 0([er fe T fe ]) C [a, 0], 

and 

P - a = (f>(T k ) - <f>{a k ) = { 0(t)df < M(r k - a k ), 



and so r k — & k > ^r. Moreover, since 4>i[a k , Tk]) Q [a, 0] it is 



2 1 oi/fj^nw ^ 



eo < e(t) = (<f>(t)) 2 + 2V(0(t)) < (<X*)) + 2F, 
for any ie [cr*, Tk]. Then 

(0(t)) 2 > e - 2K > 

for any t E [(T k ,T k ). 
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Therefore, by fl2.iq ), 

yfttjkj- V^K) = f " ^r dt = ^ f " ^ dt > (eo-2K)V3(r fc -<7 fe ) > (e -2^)V3^, 

@ k Z V e Ja k 

and then e(t) must diverge. 



Now, recalling (2.3), let > be a sufficiently small constant, and let — > i s a sequence such 
that </>(r fc ) = <f>, with u(^) < 1 - 6 as > <p. Then 0(r fc ) 2 = e(r k ) - 2V(<p) -> +oo. Let us 
suppose that <j>(rk) is unbounded above - analogously one can argue if it is only unbounded below. Up 



to subsequences, we can suppose that (f>(rk) positively diverges. Moreover, recalling (3.3), p(rh) — ► 
and then, if k is sufficiently large, 

ip^K^y-i, 

so that, since cj){rk) — 4> an d u(<f>) < 1 — 6, we have 



«(^)Wi + pW-ko. 

Then, p(rfe) < 0, that is p(t) is decreasing at r k - But we can observe that is increasing in Tk, 
ensuring u(<fi(t)) < 1 — 6, and so, in a right neighborhood of r k , 



u(<j>(t))y/l + p{t)- 1< (1 - + - 1 

that is decreasing. We conclude that the function p(t) = 2V ^t^ decreases for t > r k , until t equals 

a local maximum t k , where cb vanishes. Then there exists n > such that 2V(<j>(t)) < K(f>(t) 2 until 
t equals a local maximum t k . But this would imply e(<j)(tk)) = 0, getting a contradiction. Then is 



eventually non zero. This implies that there exists lim \<fi(t) \ so, by Proposition 3.2, it is +oo. □ 

t — >+oo " " 



We have shown so far that <f>{t) diverges, as t approaches t s = sup I. In the following, we will 
show that t s € II for almost every solution, in the sense that there exists a set of initial data, whose 
complement has zero Lebesgue measure, such that the solution is defined until a certain finite comoving 
time t s (depending on the data). 

Henceforth we will suppose (just to fix ideas) <fi(t) positively diverging and cf>(t) > Vi > i . Then, in 



the interval [£, +oo), p can be seen as a function of 4>, that satisfies by (3.4) the ODE 
(3.5) ^ = 2^3(1 + P) («(0) VT+P - l) • 

Notice that considering cf) as positively diverging is not restrictive, since it can be easily shown that p(~<j>) 
satisfies the same equation as p. 

We now state the following crucial result: 

Lemma 3.5. Except at most for a measure zero set of initial data, the function p(t) goes to zero for 
t — * sup I. 

To obtain the proof we need a result of of existence, uniqueness, and instability of solution for a 
particular kind of ordinary differential equation of first order to which the standard theory does not apply 
straightforwardly. The proof of such a result is based on a suitable estimate of the involved neighborhoods 
when applying Implicit Function Theorem. 

Theorem 3.6. Let us consider the Cauchy problem 

h(s)z(s)=f(s,z(s))+g(s), z(0) = 0. 

and (3 is a positive constant such that the following conditions hold: 

(1) h e C°([0,/3],R) such that h(0) = 0, h(s) > in ]0,/3], and /i(s) _1 is not integrable in]0,f3]; 

(2) g <E C°([0,P],R) such that g(0) = 0; 
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(3) / e C 0,1 ([0, 13} x R, R) rac/z that f(s, 0) = 0, §£(0, 0) < 0, and 

df 



(4) 3p > 5mc/z f/zaf i^(s, z) is uniformly Lipschitz continuous with respect to z in [0, 0\ x [— p, p], 



f/iaf /i > swc/z that, if\zi\, \z%\ < p, s € [0, /3], f/zen |§j(s, Zi) — ^(s, < — 22I- 
77ze«, there exists a < (3 such that the above Cauchy problem admits a unique solution z(s) in [0, a]. 
Moreover, this solution is the only function solving the differential equation h(s)z(s) — /(s, z(s)) +g{s) 
with the further property liminf s ^o \ z ( s )\ = 0. 

Proof. Let us set a < free to be determined later, and let X a be the space 

X a = {ze C°([0, a}) n C x (]0, a]) : z(0) = 0, lim Ws)i(s) = 0}. 

s->0+ 

It can be proved that X a is a Banach space, endowed with the norm ||z|| a = ||z||oo + ||ft.i||oo- 
Let also be 

y a = {AeC°([0,a]) : A(0) = 0} 
a (Banach) space endowed with the L°°-norm and let us consider the functional 
T:X a ^ y a , F(z)(s) = h(z)z(s) - f(s, z(s)). 
It is easily verified that T is a C 1 functional, with tangent map at a generic element z € X a given by 

(d^(zm)( S ) = h(s)i(s)-^( S ,z(sMs), 

where £ £ X a . Observing that g € y a , we want to find a such that the equation 

(3.6) T(z) = g 

has a unique solution z € X a . To this aim, we will exploit an Inverse Function scheme, and we will 
prove that J 7 is a local homeomorphism from a neighborhood of zq = in X a onto a neighborhood of 
T(zq) = in y a , that includes g. This will be done using neighborhoods with radius independent of a 
and this will be crucial to obtain the uniqueness result. 

In the following we review the classi c sc heme (see for instance |Jl|]) for reader's convenience. Let 
R(z) = T(z) - Tiff) - 6uF(0)[z\; then tfufy is equivalent to find z such that R(z) + AF(0)[z] = g, and 
therefore, if <£F(0) is invertible, to prove the existence of a unique fixed point of the application T on 

(3.7) T(z) = (dF(0))- 1 [g-R(0,z)}. 

We first show that T is a contraction map from the ball B(0,S) C X a in itself, provided that 6 and \\g\\oo 
are sufficiently small (independently by a). The following facts must be proven to this aim: 

(1) there exists a constant M, independent on a, such that, ~iz\,Zi € X a with ||zi|| Q , H^IU < 1, it 
is 

(3.8) ||dF(«i) - dF(^)|| < M ||zi - « 2 || Q 

(the norm on the left hand side refers to the space of linear applications from X a to y a ). 

(2) there exists a constant C, independent on a, such that 

(3.9) HdJ^or 1 !! < C 

(here the norm refers to the space of linear applications from y a to X a ). 
If the above facts hold, given z\, z 2 <E X a , then 

AT(Q)[T( Zl )-T(z 2 )] 

= R(z 2 ) - R(zi) = T{z 2 ) - T( Zl ) - dT(0)[z 2 - Zl ] 



f (dT(tz 2 + (1 - t)z x ) - dF(0)) [Z2 - zi] dt, 
Jo 
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hence, if in addition zi, z-x € B(0, 8), 
\\T(z 1 )-T{z 2 )\\ a < 

IKd^O))- 1 !! (J || {AT{tz 2 + (1 - t) Zl ) - dF(0)) || d?j \\ Z1 - z 2 ||oo < 2MCS\\ Zl - z 2 \\ a . 
and therefore: 

• T is a contraction, taking 8 such that K := 2MC8 < 1; 

. since ||r(«)|| Q < \\T{z) - T(0)|| Q + ||T(Q)|| a < tf||z|| a + || (cLF(O))- 1 [<?] || a < K\\z\\ a + 
Cllfflloo, then T maps B(0,8) in itself, provided that H^H^ < C _1 (l - Jf)<5. 
Observe that the first of these two facts fixes the value of 8, whilst the last inequality holds choosing a - 
free so far - small enough. This is one of the reasons why the constants M and C must be inde pend ent 



on a. Then T admits a unique fixed point on B(0, 8), which is the solution to our problem (^6) on 
the interval [0,a]. In other words, the function T is a local homeomorphism from B(0,6) C X a to 

s(o,c- 1 (i--K')«)cy a . 

To see that the solution is globally unique on X a , let us argue as follows. Suppose z £ X a \ B(0, 8) 
(z z) solves the problem, and let a\ < a sufficiently small such that ||z|| ai < 8. T hen, observing 



l5lli°°([o,a;i]) ^ llslli^ffOja]) ^ C -1 (l — K)S, and recalling that estimates (3.8)— (3.9) do not depend 
on a, one can argue as before to find that 2|ro ai i is the unique element of B(0, S) C X ax mapped into 
Plfo.ai] G -B(0,C _1 (1 — K)S) C J^ai- But, of course, ||z|| Ql < ||z|| Q < 8, so z and z coincide on 
[0, a±], and therefore on all [0,a]. 

Therefore, to complete the proof of the existence and uniqueness for the given Cauchy problem, just 



(3.8 )— (3.9) are to be proven. The first equation is a consequence of local uniform Lipschitz continuity of 



/ )2 . The second one needs some more care: taken A s y a , we must consider the Cauchy problem 
(3.10) Hs)i(s) =£(s)£(s) + \(s), £(0)=0, 

where £(s) := j £(s, 0), that without loss of generality we can suppose negative, Vs 6 [0, 0]. First, it is 



easily seen that (3.1C) admits the unique solution £ £ X a , 



Then IKd^O))- 1 ^]^ = HfHa < (1 +4)||elU + ||A|U, where i x = ||^|Uo= [0lffl . Moreover, called 
^0 = sup[ ^ £(s) < 0, it is easily seen that ||£||oo < — ^HAHoo, so it suffices to choose C = 1 — , 



and (3.9) is proven. 

To prove last claim of the Theorem, let us suppose that w(s) is a function defined in [0, a] such that 
h(s)w(s) — f(s,w(s)) + g(s), and that Sk is an infinitesimal and monotonically decreasing sequence 
such that w(sk) — * as k — * 00. We want to prove that w = z, and therefore, it will suffice to show that 

lim s ^ w(s) = 0. 

First of all, observe that from the hypotheses, the equation f(s,z) + g(s) — defines a continuous 
function £(s) : [0, 8] — > E, such that C(0) = 0. In particular, since ||(0, 0) < 0, 3p > such that, in 
the rectangle [0, 8] x [— p, p], it must be f(s, z) + g(s) < (resp.: > 0) if z > ((s) (resp.:z < C( s ))- 

Let us now argue by contradiction, supposing the existence of an infinitesimal sequence <7fc, that can 
be chosen with the property < Sfc, such that | w(tTfc ) | > 8 for some given constant 9. Now, up to 
taking a smaller constant 8, then \G(s)\ < |, Vs £ [0, 8]. Then, for k sufficiently large, |iw(<7fc)| > 9 > 
I > sup[ 5 ] |C(s)|- Therefore, Vs < at, \w(s)\ > \w(ak)\, which is a contradiction since w(sk) — * 0. 
Then lim s ^ w(s) = 0, and the proof is complete. 

□ 



Proof of Lemma 3.5. Than ks to Proposition 3.4 we can carry on the proof by studying qualitatively the 
solutions of the ODE (3.5). We will be interested in those solutions which can be indefinitely prolonged 
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on the right. With the variable change y = yT 

dy 



p, (3.5) becomes 

= Vi{y 2 - 1) (u(4>)y - 1) , 



(3.11) 

u<; 

where we recall that u{4>) is give n by ( |2.2| ). 

Let us consider solutions of ( [3.1 1| ) defined in [<po, +oo). If liminf^ +00 \y(<f))u(<f>) — 1| > then, 
necessarily, u((f>)y(<f)) — 1 < (and y((f>) decreases, and goes to 1), otherwise the solution would not be 
defined in a neighborhood of +oo. In short, if y((f>) is a solution defined in a neighborhood of +oo, the 
following behaviors are allowed: 

(1) either y((f>) is eventually weakly decreasing, and tends to 1 (so that p tends to 0), 

(2) orliminf^+oo \y(4>)u{<j>) - 1| = 0. 
With the variable and functions changes 

s = eT^^ , z = uy — 1, 



equation (3.11) takes the form 
(3.12) (su)z(s) 



-(z + l) 2 z + (u 2 + s u)z + s u. 



But it can be easily seen, using (2.4), and the identity 

su(s) = -^u'(<f>), 



that ( 3.12 ) satisfie s the assumption of Theorem p.q , and so there exists a unique solution for the Cauchy 
problem give n by ( 3.12 ) with the initial condition z(0) = 0, that is furthermore the only possible solution 
of the ODE ( 3.12 ) with liminf s _,o \ z ( s )\ — 0, and this fact results in a unique solution satisfying case 
(^) above, whereas all other situations lead to case ([!]). □ 

Now, we are ready to prove the main result: 

Theorem 3.7. IfV{4>) € € then, except at most for a measure zero set of initial data, there exists 
t s G R such that the scalar field solution becomes singular at t = t s , that is liiru^- e(t) = +oo, and 

linL^- a(t) = 0. Moreover lim.^- <j>{t) = +oo. 

Proof. We divide the proof into two parts, depending on the value of the initial energy. 
*) case e > 2V^. 



From ( |2.7D and ( |2.9a| )-( |2.9b[ ), it easily follows that 
(3.13) e(t) = 2\/3e 3/2 



1 



1 + P(t) 



But Lemma 3.5 ensures that p (defined in ( |3.3| )) goes to zero for almost every choice of initial data. For 
this choice, e(t) > e 3 / 2 in a left neighborhood of supl. Then t s < +oo and lim t _ >t - e(t) — +oo easily 
follows, using comparison theorems in ODE. To prove that lim f _ i . t - <fr(t) = +oo, it suffices to observe 
that 

e(t) 



Pit) + 1 



1. 



m 2 

0, we first prove that 
lim 

/o 

Indeed, recalling that pit) is eventually bounded, 3i > such that e(t) — cf) 2 (t)(l 
for some suitable constant n > 0, so y/e(t) > K(f>(t), Vt > t, and then 

ft r t 



Finally, to show that lim t ^ t - a(t) 
(3.14) 



Jo 



p(t)) > K ^(tf 



y/edr > k 
t Jt 



dr = K(4>(t) - (f>(t)). 
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The righthand side above diverges, so that (3.14) is proved; but, using (2.9a), we get 



log 



ait) 
o(0) 



/ -dr= — [ y/edr, 
Jo a V3 Jo 



from which lirru^- a(t) = follows, since by the righthand side is negatively divergent. 
ii) case < e < 2\4. 

Suppose eo = 0. If the initial data are not at the zero level of V (which consists only of a finite 
number of points because of assumption (Al)), thanks to Lemma B.l in Appendix [B], where a result of 
local existence and uniqueness for the solutions of equation ( 2.9b| ) satisfying e(0) = 0, but such that 
e(t) > for t > 0, we can assume eo > without loss of generality. 

Now suppose that < eo < e(t) < 2V* Vt £ I (otherwise the previously obtained results would 
apply). We shall show that the solutions for which this holds are non-generic. 

To show this fact we observe that in this case (j>(t) £ [a, (3} for any t > and, since is bounded, 
sup I = +oo. Then Lemma A. 1 in Appendix ^ applies and we have 



lim 4>{i) = 0, lim V (<£(*)) = 0. 

t — >+oo t — >+oo 

Moreover by assumption (Al) there exists 0* critical point of V such that lim^ +oc <f>(t) = tf>*, (with 
critical value V(</>») = lim t ^ +oc e(i) G]0, 2T4].) Whenever is a (nondegenerate) maximum point 
we can study the linearization of the first order system equivalent to ( 2.9b ) in a neighborhood of the 
equilibrium point (0, 0*) (such that V{4>*) > 0) , as done in JTlh, obtaining that there are only two 
solutions cf>i and <f>2 of (2.9b) such that lim <fii(t) — (f>*, lim <fii(t) = 0, i = 1,2. Moreover there 

t — >+oo t — >+oo 

exists t* > such that <fii (t) > for all t > t» and (j>2 (t) < <j>* for all t > £*. It can be shown Jill ] that, 
if 0* is a minimum point and <j) starts with initial data close to (0, </>*), then <\> moves far away from 0*. 
Therefore, under the assumptions made, for almost every choice of the initial data, the function <f>(t) must 
be such that its evolution cannot be contained in the bounded closed interval [a, /3], namely, e{t\) > 214 
for some t\ > 0, which allows us to apply the theory we already know to show that the singularity forms 
for almost every choice of the initial data. □ 



Remark 3.8. Since the potential can be negative the energy e, could in principle vanish "dynamically" 
(namely <f> could be a solution of (/>(t) 2 = —2V(4>(t)), Vi). However, these functions are not solutions of 
Einstein field equations ( 2.5a )-( 2.5b ), since e(i) = implies a(t) = ciq and then cj)(t) — 0, V((j)(t)) = 
V(<po) = 0. Thus, one is left with a set of constant solutions with zero ener gy of the fo rm ( ap, <f>o), with 
ao positive constant and V(<f>o) = 0. By the results obtained in Lemma B.l in appendix B.l we therefore 
see that at the "boundary of the weak energy condition", local uniqueness of the field equations may be 
violated if V'((j>o) ^ and V((j>) = 0. However the set of the initial data for the expanding equations 
intersecting such solutions in a finite time has zero measure (note that the points </>o such that V(4>q) = 
and V'((j)o) ^ are isolated). 

3.1. Examples. 

Example 3.9. The above results hold for all potentials with polynomial leading term at infinity (i.e. 
\ 2 <fi 2n ). For instance, for a quartic potential V{<p) = —^m 2 <fi 2 + A 2 </> 4 , with A,m/0, the function u(<f>) 
goes as — %- for <f> — > +oo, and all conditions listed above hold. The behavior for a particular example 
from this class of potential is given in Figure |1| 

Example 3.10. In the case of exponential potentials, the results hold for asymptotic behaviors with lead- 
ing term at infinity of the form Vbe 2v/ ^ A I^I with A < 1. Consider for instance V(<f>) = Vo cosh(2v / 3A0), 
where Vo, A > 0: the quantity u(<p) goes like A, and so (2.3) is verified if A < 1. See a particular situation 
from this class represented in Figure [| 
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Figure 1. Behavior of the function <fi(t),e(t), and a(t) with potential given by 
V(<fi) = 1 — 4> 2 + </> 4 . The initial conditions are 4>q = —0.6, ct>o = 0, ao = 1. 
The time of collapse t s approximately equals t s = 2.1. 




FIGURE 2. Same as Figure 111, but now the potential is given by V(<fi) = 
cosh(-\/60)/lOO. The initial conditions are cf>o = 0.1, <po = 0, ao = 1. The time 



of collapse t s approximately equals t s = 18.6. 



Example 3.11. Decaying exponential potentials can be considered as well. For instance for V(<p) — 
(1 — e~ Q V'* 2 +T 2 ) 2 the function u(<f>) behaves like e _Q '^', so it goes to zero at ±oo. Figure ^| represent 
a situation where some corrections terms have been added in order to obtain a potential with more than 
one critical points. Even choosing initial data such that eo = 0, the solution diverges in a finite amount 
of time, as showed in general in the second part of the proof of Theorem 3/7. 



4. Gravitational collapse models 

In what follows, we construct models of collapsing objects composed by homogeneous scalar fields. 
To achieve this goal we must match a collapsing solution (considered as the interior solution in mat- 
ter) with an exterior spacetime. The natural choice for the exterior is the so-called generalized Vaidya 
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eft) 




FIGURE 3. Same as Figure where now V{4>) = 2(1 - e -vWT)2 _ i e -(*+4) 2 _ 
l e -(i-4) 2 _|_ ^ -p ne initi al conditions are (f> = 0, </> = \/— 2V(0), a = 1. The time 
of collapse t s approximately equals t s — 1.8. 



solution 

(4.1) d.s c 2 xt = -(i- 2M( ^' y) j d[/ 2 -2drdt/ + y 2 dr! 2 , 

where M is an arbitrary (positive) function (we refer to [ |27| ] for a detailed physical discussion of this 
spacetime, which is essentially the spacetime generated by a radiating object). The matching is performed 
along a hypersurface S which, in terms of a spherical system of coordinates for the interior, has the simple 
form r = r b = const. The Israel junction conditions at the matching hypersurface read as follows [^]: 

(4.2) M(U(t), Y(t)) = ^r 2 b a{t)a(t) 2 , 

(4.3) W {U(t) > Y(t)) = \ r " {h{t)2 + 2a W a (*))> 
where the functions (Y(t), U(t)) satisfy 

(4.4) Y(t)=r b a(t), ^(t) = ' 



dt 1 + a(t)r b ' 



The two equations ( 4.2 >— ( )4.3| ) are equivalent to require that Misner-Sharp mass M is continuous and 
Qfij- = on the junction hypersurface S (using ( |4.2[ )-( |44[ ) together with ( 2.5a| )— (p~^), it can be seen that 



in this way also the equation of motion for the scalar field remains smooth on the matching hypersurface). 

The endstate of the collapse of these "homogeneous scalar field stars" is analyzed in the following 
theorem. 

Theorem 4.1. Except at most for a measure zero set of initial data, a scalar field "star" collapses to a 
black hole. 



Proof. It is easy to check that the equation of the apparent horizon for the metric (2.1) is given by 



r a{t) = 1. If a is bounded, one can choose the junction surface r — r b sufficiently small such that 
a 2 (i) < -pr, V(t, r) 6 [0, t s ] x [0, rb], and so (l — i s bounded aw ay fr om zero near the singularity. 



As a consequence, one can find in the exterior portion of the spacetime (4.1), null radial geodesies which 
meet the singularity in the past, and therefore the singularity is naked; otherwise, if d 2 is unbounded, the 



trapped region forms and the collapse ends into a black hole [Bfl. Now, we have proved in Lemma 3.5 the 
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existence of p^ = lim^oo p{4>) (which is also equal to lkcu^- pit)), and that, actually, poa vanishes 



except for a zero-measured set of initial data. On the other hand, using ( 2.5b ), we get 
,4,) * = --( 2 TT 

and therefore a < — — eventually holds in a left neighborhood of t s , say [to, t s [, where d is decreasing. 
It follows that \/t G]ioi ts [, o,(t) < a(to), which is negative by hypotheses. Then 

d(t)-d(i )= / d(r)dr<-/ ^Vdr<-d(i )/ 44 dr > 
Jt Jt a(r) J to a(r) 

that diverges to — oo. Then a(t) is unbounded. □ 

5. Naked singularities as examples of non-generic situations 
As we have seen, the formation of naked singularities in the above discussed models is generically 



forbidden by theorem 4. 1 . Thus, the theorem does not exclude the existence of initial data which give 
rise to naked singularities for the given potential, it only assures that such data are non-generic. Actually, 
there are examples of such solutions containing naked singularities in the literature, and it is our aim here 
to show how they fit in this scenario. 

To discuss this issue, it is convenient to reformulate the dynamics in the following form [Bp. Since 



a(t) is strictly decreasing in the collapsing case, we re-write equations (2.7), (2.9c) and (2.9b) using a 
instead of t as the independent variable. Setting 



(5.1) e(o) = 3 

we get: 

(5.2a) d = —if) (a), 



2 



1 / 



(5 - 2b) [fa) =^{ 1 -%(a) 

(5.2c) V=(m 2 (l + la^ 



2 ip(a) 

In the above equations, V(<p) plays of course the same role of "equation of state" as before. However, one 
can try to prescribe a part of the dynamics and then search if there is a choice of the potential and a set 
of data which assure that this dynamics is actually an admissible one for the system at hand (of course, 
in this way, only the "on shell" value of the potential can be reconstructed). In particular, a somewhat 
natural approach can be that of prescribing the energy (e(a) or equivalently "J (a)). 

In jl2[], the case in which the energy density behaves like aT v with v positive constant is considered. 
It is then shown that for v < 2 the system forms a naked singularity, and the existence of an interval of 
values of v is interpreted as a generic violation of cosmic censorship (that is anyway restored when loop 
quantum gravity modifications are taken into account, see [|l3|]). However, due to the equation ( |5.2c| ) 
above, the on shell functional form of the potential for this example can be calculated to be V((j>) = 
Vbe 2v/ ^ A l^!_with A = \jvj§. This potential belongs to the class which has been treated in our previous 



example p.lOt and therefore we know that a generic choice of data must give rise to a blackhole (for this 



class the function u(<ft) = V /{2\f?>V) behaves like a strictly positive constant A < 1). However, as 



we know from the proof of Lemma [3.5[ , a non generic situation can occur, when lim t ^ ( - uy/1 + p = 1 
(where p = 2V/(p 2 ); which means that lim^oo p(<f>) = 4$ — 1 = 6/v — 1. Evaluating the function p 
for the solution under study, we get exactly p — Q/v — 1. Therefore, the case discussed in [|l2| [jj]] is 
precisely the non-generic one. The constraint on the data which implies non-genericity comes, of course, 
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e(t) 



l-2M/Y(U(t),Y(l» 




FIGURE 4. An example of naked singularity for V(<fi) — O.Ole* for 4> > 0. The initial 
conditions are <j>o = 0, </>o = 550 -1 / 2 , ao — 1. The matching with external solution is 
performed at = 100. The time of collapse t s approximately equals t s = 38.8. 



from the fact that imposing a certain behavior on the function ip(a) selects a two-dimensional set in the 
space of the data, characterized by the algebraic constraint </>(0) 2 + 2V((f>(0)) = 3(tp(a(0)) / a(0)) 2 . 
To complete the analysis, we observe that, although the evolution of non generic data is not predicted 



by Theorem 3.7, we can easily obtain information also in this case. To this aim, observe that, since p 



is finite, equation (3.13) can be used to prove singularity formation in a finite amount of comoving time 
with an argument similar to that used in the proof of Theorem [3.7[ Moreover, the metric function a(t) 
behaves asymptotically as a (t s — t) 1 /^ 3 * ); then, if A < \/3/3, the function a(t) is bounded near the 
singularity, confirming naked singularity formation (this choice of A corresponds exactly to the choice of 
v < 2 in ml). 

Figure ^sketches an example of this kind: the initial data are chosen in the non generical way such 
that p goes like a positive constant as t — > t s . The singularity forms (i.e. e(t) diverges in a finite time), but 
choosing rt, not too large, the function 1 — 2M/Y on the junction hypersurface r — r& remains strictly 
positive along the evolution up to the singularity. 

The above example suggests that hypotheses on the asymptotic behavior of the function u((j>) can be 
added, to obtain singularity formation also in the non generic case: indeed, if the potential V((f>) satisfies 
the condition 

(5.3) liminf | | > 0, 

(p — »±oo 

i.e. u(4>) is definitely bounded away from zero, then p is finite also in the non generic case, and equation 



(3.13) in the proof of Theorem 3.7 can be used to prove the formation of the singularity in a finite amount 



of time. To analyze the nature of this singularity observe that, if 
(5.4) liminf \u(<j))\ > 

d>— >±oo 



the argument in the proof of Theorem 4. 1 applies, and therefore the collapse ends into a black hole. On 
the other side, if 

(5.5) limsup|u(0)| < ^, 

d> — *±oo 
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then a is eventually positive, and this implies that a is bounded, so that the formation of the apparent 
horizon is forbidden. As a consequence, the collapse ends into a naked singularity. 

6. Discussion and Conclusions 

We have discussed here the qualitative behavior of the solutions of the Einstein field equations with 
homogeneous scalar fields sources in dependence of the choice of the self-interacting potential. It turns 
out that whenever a potential satisfies a certain set of general conditions singularity formation occurs for 
almost every choice of initial data. Matching these singular solutions with a Vaidya "radiating star" exte- 
rior we obtained models of gravitational collapse which can be viewed as the scalar-field generalization 
of the Oppenheimer-Snyder collapse model, in which a dust homogeneous universe is matched with a 
Schwarzschild solution (the Schwarzschild solution can actually be seen as a special case of Vaidya). 

The Oppenheimer-Snyder model, as is well known, describes the formation of a covered singularity, 
i.e. a blackhole; it is actually the first and simplest model of blackhole formation ever discovered. The 
same occurs here: indeed we show that homogeneous scalar field collapse generically forms a blackhole. 
The examples of naked singularities which were found in recent papers [j~2|], turn out to correspond 
to very special cases which, mathematically, are not generic. Therefore, our results here support the 
weak cosmic censorship conjecture, showing that naked singularities are not generic in homogeneous, 
self-interacting scalar field collapse. 

Non-genericity was already well known for non self-interacting (i.e. V((f>) = 0) spherically symmetric 
scalar fields. Whether the results obtained here can actually be shown to hold also in the much more 
difficult case of both inhomogeneous and self-interacting scalar fields remains an open problem. 

Appendix A. Some properties of global solutions 

We give here the proof that velocity and acceleration of solutions which extend indefinitely in the 
future must vanish asymptotically, if the values of the derivatives of the potential remain bounded on the 
corresponding flow. 



Lemma A.l. Let <p(t) a finite-energy solution of (2.9b), with V G £ Then: 

1) 4> can be extended for all t > 0, 

2) ifV is bounded on the flow, then lim^+oo <j>(t) = 0, 

3) if also V" is bounded on the flow then lim t ^ +00 (p(t) = 0. 

Proof. Since e is bounded and V is bounded from below, <j) is bounded as well, and thus <fr can be extended 
for all t (actually the existence of a lower bound for V is the unique condition among Al)-A3) used here 
and below). 

Using ( 2.10| ) we have 

<j){t) 2 dt = I — y ' dt < +oo, 



Jo 

so there exists a sequence tk — > +oo such that 4>{tk) — > 0. By contradiction, suppose that 3p > 0, 
and a sequence s& — ► tk - that can be taken such that Sk > tk - with |^(sfc)| > p. Let k such that 
|</>(ifc)| < f , Vfc > fc, and let rfe, sequences such that tk < Tk < Oft < and 

I0fa)l = f, \i>(a k )\=p, P - <0(hfe,o*])| <p. 



Note that e bounded implies <j> bounded, while by assumption, V'(4>(t)) is bounded. Therefore by (2.8b) 
there exists M > such that 

\4>{t)\ < M Vt. 

Therefore 

P 



2 



< \4>{a k ) - 4>{r k )\ < / \4>(t)\ dt < M(a k - r fe ) 
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that is <7fc — Tfc > 2Xf, and therefore 



r+°° r ak / r,\ 2 ^,3 

/ m 2 dt>j2 ^) 2 di>E(^-^)-(f) >E^ 
j ° k jTk k k 



that diverges. This is a contradiction, and therefore it must be lim t ^ +00 cb(t) = 0. 

To prove that the acceleration also vanish, let us first observe that 3tk — > +oo such that V (4>(tk)) — > 
- otherwise, there would exists k > such that |V'(0(t))| > k definitely, which would imply, in view of 



(2.8b), that |0(i)| > k/2, that is absurd since lim t ^ +oc <p{t) = 0. 

Then, let us suppose by contradiction the existence of a constant p > 0, and a sequence Sfc — > +oo 
such that ifc < s& and | V (0(sfc))| > p. Therefore, one can choose Ofc, Tfc sequences such that tk < r& < 
(Tfc < Sfe, and 

m#7*))| = f ■ \V'(<f>(<Tk))\ =P, | < |V'(0(t))| < P, Vi G [Tfc, (Tfc]. 

Then, since by assumption V"(cf>(t)) is bounded, there exists a constant L > such that 

f = 1^ ~ V (0(Tfc))| < jT |F"(0(i))| |0(*)|di < L(<7 fc - Tfc). 



But for sufficiently large k let us observe that ( 2.8b ) implies |0(i)| > p/4, Vt G [rfc, Ofc], and therefore 
|0(cfc) — ^(^fe)! = | J^* 4>(t)dt\ > that is a contradiction since 0^0. Thus lim f ^ +oc V'(0(t)) = 



vanishes, and equation ( 2.8b ) implies that also vanishes in the same limit. 

□ 

Appendix B. Local existence/uniqueness of solutions with initial zero-energy 

Lemma B.l. Let 0q, vq such that Vq + 2V{4>o) — 0. Then, there exists t* > such that the Cauchy 
problem 



4>(t) = -v'W)) + ^3(0(i) 2 + 2^(0(i)))0>), 

(B.l) \ 0(0) = 00, 

0o = vq, 

has a unique solution <fi(t) defined in [0, t„] vw'f/z f/ze property 

(B.2) e(t) = s( [ <j>(s) 2 ds) , Vt G]0,fc*]. 



vo / 

Moreover if (0o,mi ^o,m) ~~ * (0o, w o)> (^o.m ) 2 + 2V(0o ;m ) = ant/ m is f/ze solution of (BT) wz'f/z 
initial data (0o,m,^o,m) satisfying condition ( B.2 ), z'f z'i m — > wz'f/z respect to the C 2 -norm in the 
interval [0, t»]. 

Proof. Let us consider the "penalized" problem 

[ 0(<) = -F'(0(*)) + V^W" + 2y (^)) + 5) 
(B.3) j 0(O ) = 0,,, 

[ 00 = uo, 

that has a unique local solution 0„. If n is not defined Vt > 0, let I„ be the set 

I„ = {t G R : |0„(s)| < |0 O | + 1, |0„(s)| < |vq| + 1, Vs > i}. 
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Of course, I„ ^ and, called t n = sup In, if t n is finite, then \cb n {t n )\ = |«o| + 1, or \(j) n {tn)\ = |0o| + l- 
Now assume \4> n (t n )\ = \4>o\ + 1- Then 

1 = |0n(*n)-<fo| < / \Ms)\ds<(\v \+l)t n . 

Jo 

Analogously if \(f> n (t n )\ = \ v o\ + 1 we nave 

1 = \<j>n(tn) ~ V \ < / |^n(s)|ds. 

Jo 



Since |<75 n (£)| < |</>o| + 1 an d I0n( s )l < \ v o\ + 1 f° r all t £ [0, i„], and qb n solves (B.3), we see that there 
exists C independent of n such that \4>{t)\ < C for all t G [0, t n ] ■ Therefore in this second case we obtain 
1 < Ct n . 

Then t* := inf„ t n > (we set t* = 1 if <„ = +oo Vn). Moreover \cf> n \ is uniformly bounded in 



[0, then up to subsequences, there exists a C 1 function <r3(t), solution of (B.l ), such that <p n — > and 
(73„ — > <j) uniformly on [0, t*]. 

Now, consider e n (t) := <j) n (t) 2 + 2V{<j> n {t)) + ^x- We have 

(B.4) e n {t) = 2VS^(t)Mt) 2 . 

Then e is not decreasing, while e(0) = ^. Then is uniformly bounded away from zero and therefore by 



(B.4), dividing by ^f&H and integrating gives y/e n (t) = + V% J Q <Pn(s) 2 ds. Therefore passing to the 

limit in n we obtain e(t) = <ft(t) 2 + 2V(4>{t)) = 3 ( f* 4>{s) 2 ds^j for all t € [0, t»] obtaining the proof 
of the existence of a solution. 

The uniqueness of such a solution can be obtained by a contradiction argument. Assuming <f> and tp 



solutions, and called 9 = tp — ip, one can obtain, using (B.l ), the estimate 



\0(t)\ <Kt f \6(s)\ds + K 2 f \e(s)\ds 
Jo Jo 



for suitable constants K\, Setting p(t) = \0(t)\ + \0(t)\, and observing that p(0) = 0, it is not hard 
to get the estimate p(t) < (K\ + K 2 + 1) L p(s) ds, and then p = from Gronwall's inequality. 

Finally using Gronwall's Lemma as above we obtain also the continuity with respect to the initial 
data. □ 

Remark B.2. Reversing time direction in the above discussed problem ( |B.1[ ) yields a results of genericity 
for expanding solutions such that the energy e(t) vanishes at some finite time T, 
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